Quark Spectral Function and Deconfinement at Nonzero Temperature 
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The maximum entropy method is used to compute the quark spectral function at nonzero temper- 
ature. We solve the gap equation of quantum chromodynamics (QCD) self-consistently, employing 
a rainbow kernel which phenomenologically models results from Dyson- Schwinger equations (DSE) 
and lattice QCD. We use the criterion of positivity restoration of the spectral function as a signal 
for deconfinement. Our calculation indicates that the critical temperature of deconfinement Td is 
slightly smaller than the one of chiral symmetry restoration T c : Td ~ 94%T C in the chiral limit, and 
Td ~ 96%T C with physical light quark masses. Since these deviations are within the systematic error 
of our approach, it is reasonable to conclude that chiral symmetry restoration and deconfinement 
coincide at zero chemical potential. 



I. INTRODUCTION 

Heavy-ion experiments at the Relativistic Heavy-Ion 
Collider (RHIC) and the Large Hadron Collider (LHC) 
are focusing on charting the phase diagram of hot and 
dense nuclear matter. The quark-gluon plasma (QGP), 
a primordial state of matter in the early Universe, where 
chiral symmetry is restored and quarks and gluons are 
deconfincd, has been re-created in the extremely hot en- 
vironment of a heavy-ion collision. With the expansion of 
the fireball, nuclear matter cools down and dilutes. The 
low-temperature, low-density phase of nuclear matter is 
characterized by confinement and dynamical chiral sym- 
metry breaking (DCSB). It is a central goal of modern 
theoretical physics to understand the properties of, and 
the transitions between, these phases. The chiral and 
deconfinement phase transitions as well as their inter- 
play are especially interesting. In general, the chiral con- 
densate (or, equivalently, the dynamical quark mass) is 
adopted as an order parameter for the chiral phase tran- 
sition. The order of this transition may depend on the 
number of quark flavors, the values of the quark masses, 
and whether the U(1)a anomaly of QCD is effectively 
restored. For nonzero quark masses, the existence of a 
critical end point (where the transition turns from be- 
ing first order at low temperatures and high densities to 
being crossover at high temperatures and low densities) 
has been suggested but, even if it exists at all, its precise 
location is still highly debated. 

Concerning the deconfinement phase transition, the 
situation is even more complicated because confinement 
has been a mystery since the inception of the Standard 
Model. The notion of confinement is easily understood 
from the linearly rising potential between infinitely heavy 
uarks [l],[2|, which has also been studied by lattice QCD 
However, this is no longer true for light quarks 
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because of strong pair-creation and -annihilation effects 
|J|. In the pure-gauge limit realized for infinitely heavy 
quarks, the center Z(3) symmetry of the color gauge 
group SU(3) is preserved in the confining phase, while it 
is spontaneously broken in the deconfined phase. Here, 
the Polyakov loop (or the thermal Wilson line) 0, Q 
is the corresponding order parameter. Equivalently, the 
dual quark condensate [7|-|9[ was proposed as an order pa- 
rameter, which makes it possible to study the interplay 
between confinement and DCSB. However, their validity 
as order parameters for phase transitions in light-quark 
systems remains unclear. 

Besides these order parameters, confinement can be re- 
lated to the analytic properties of QCD Schwinger func- 
tions [10l4l2| . The axiom of reflection positivity requires 
that the propagator must have a positive definite Kallen- 
Lehmann spectral representation for asymptotic (or de- 
confincd) quarks. In other words, if the quark propaga- 
tor can be decomposed in terms of complete eigenstates 
of the Hamiltonian, each of which should have positive 
probability, quarks can propagate as asymptotic states; 
otherwise, quarks have to be somehow confined. It can 
be shown that pairs of complex conjugate poles of the full 
quark propagator lead to a nonpositive definite Kallen- 
Lehmann spectral representation. Therefore, in Refs. 
[llj 13 ; the existence of such pairs of complex conjugate 
poles was considered as a criterion for confinement. How- 
ever, this is only a sufficient condition. The reason is that 
the violation of reflection positivity can also be realised 
by propagators with real poles, namely if they have a 
negative residue which also leads to the Kallen-Lehmann 
spectral representation being not positive definite [la ]. 
On the other hand, the positivity of the quark spectral 
function is a necessary and sufficient condition for quark 
deconfinement, no matter whether the singularities of the 
quark propagator are located on or off the real axis. In 
short, by considering the quark spectral function directly, 
one is able to distinguish confined phases (where reflec- 
tion positivity is violated) from deconfined ones (where 
reflection positivity holds). 

In this work, we use the maximum entropy method 



(MEM) [16T - KL9| to explicitly compute quark spectral 
functions from the self-consistent numerical solution of 
the QCD gap equation. We employ a rainbow kernel 
J20L l2l| which phenomenologically models recent results 
from DSE [H, [H] and lattice QCD [Hl^. We define 
a deconfinement temperature as the temperature above 
which the positivity of the quark spectral function is re- 
stored. This work is a continuation of Ref. [27| where 
the quark spectral functions were studied in the region 
above T c . This paper is organized as follows. In Sec. 
III1 we present the QCD gap equation and discuss the 
ansatz employed for its solution. In Sec. IIII1 we derive 
the relation between the quark spectral function and the 
solution of the gap equation. Here, we also define the 
order parameter which signals the positivity of the spec- 
tral function. In Sec. IIV1 we briefly outline the MEM 
and its extension for nonpositive definitive spectral func- 
tions. Section [V] reports our numerical results. Finally, 
we conclude with a summary and some remarks. 



II. QCD GAP EQUATION 

At nonzero temperature, the QCD gap equation is 

Siiujn^y 1 = tf ■ p + i-fiUJn +m + T,(iuj n ,p), (1) 

E(ico n ,p) = — Y^ / 7^3 9 2 D^(k,n n i) 
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xj fl S{iuj l ,q)T u (q,uj l ,p,uj n ) , (2) 

where ui n = (2ra + l)7rT is the fcrmionic Matsubara fre- 
quency, k = p — q, Q n i = Lu n — loi, D^ is the dressed 
gluon propagator, and r„ is the dressed quark-gluon ver- 
tex. The solution of the gap equation can be expressed 
as 

Siiu^p)- 1 = ij-pA(ujl,p 2 ) 

+ i li u Jn C{ojl,p 2 ) + B{ujI,p 2 ) (3) 

or, equivalently, 

S(iuj n ,p) = - i^ ■ p a A (u? n ,p 2 ) 

-ij4U n <Jc(ul,p 2 ) + <JB{ul,p 2 ), (4) 

where A,B,C, and <ja,b.c arc scalar functions. The 
dynamical quark mass is defined as M(yi 2 ll p 2 ') = 
B(w^,p 2 )/A(oj^,p 2 ), which is independent of the renor- 
malization point. In the chiral limit, the chiral conden- 
sate is defined as 
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However, because of an ultraviolet divergence the inte- 
gral in the above equation is not well defined at nonzero 
current quark mass. Conveniently, M :— M(ujQ,p 2 = 0) 
can be used as the order parameter for the chiral phase 
transition, which is equivalent to the chiral condensate. 



The gap equation is closed by specifying the vertex and 
the gluon propagator. Here, we use the rainbow trun- 
cation, i.e., the leading term in a symmetry-preserving 
scheme 1281 : 



g 2 D tlu (k, Ct n i)T v (q, ui, p, u n ) 

= [P^D T (k 2 7 il 2 nl ) + P^D L (k 2 7 il 2 nl )}^ 
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where P^ L are transverse and longitudinal projection 
operators, respectively, 
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with T(s) = [1 — cxp(— s/^m 2 )]/ s, r — e 2 — 1, m t = 
0.5 GeV, j m = 12/25, and Aq^ 4 = 0.234 GeV. For 
pseudoscalar and vector mesons with masses < 1 GeV, 
this interaction provides a uniformly good description 
of their vacuum prop erties when aD = (0.8 GeV) 3 and 
a G [0.4, 0.6] GeV 0,[2l|, which means that there is only 
one free parameter in the model. The physical masses of 
the light quarks are vrv d = 3.4 MeV at our renormaliza- 
tion point £ = 19 GeV. Generalizing to T ^ 0, we have 
followed perturbation theory and included a Debye-like 
mass in the longitudinal part of the gluon propagator: 
m 2 = (16/5)T 2 [for details, see Ref. [Ig]. 



III. SPECTRAL REPRESENTATION 

The dressed quark propagator is related to the retarded 
real-time propagator by analytic continuation, 

S R {oj,p) = S{iui n ,p)\ iUn ^ u+ie . (11) 

From the spectral reprentation of S R (uj,p), i.e., 

p{u,p) = -2QS R (u,,f), (12) 

one immediately obtains 

S(iw n ,p)= — - -. (13) 

J_ oc 2tt iw n - u/ 

According to Eq. (j4]), the spectral function can be de- 
composed as 

P(u,p) = -n -pp v (uj,p 2 ) 

+ J 4 UPe(w,p 2 ) + p s (uJ,p 2 ). (14) 



As a consequence of the anti-commutation relation, the 
spectral function fulfills the following sum rule, 



f + °° du , ... 
/ ^-p(w,P)74 = 1 
J-oo ^tt 



(15) 



Then one can define the spectral function po(to,p 2 ) := 
cup e (oj,p 2 ), which is nonnegative and can be treated as 
a probability distribution for deconfined quarks. Note 
that po(oj,p 2 ) can be easily related to the dressed quark 
propagator by 



So(col,p 2 ) = iuj n< T C (ujl,p 2 ) , 

<- +00 dm' po(LU / ,p 2 ) 

, 27T U 1 — iuJr, 



(16) 



where oc is the scalar function in Eq. (QJ, or to the 
imaginary-time quark propagator 



D (T,p 2 ) =Tj2e- iUnT So(u 2 n ,p 2 ), 



du e (l/2~TTMT 
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p (w,p 2 ). (17) 



The above equations connect the quark spectral func- 
tion which we consider to the numerical solution of the 
gap equation. As a signal for positivity violation (or 
restoration) of the spectral function, one can define an 
"order" parameter as 
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(18) 
(19) 
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where p(ui) simply denotes po{ui,p 2 = 0). It is apparent 
that Z p = Z p = 1 and L p = for a positive definite spec- 
tral function, otherwise Z p > Z p = 1 and L p > 0. The 
critical temperature Td of the dcconfincmcnt transition 
is defined as the lowest temperature where L p = 0. 



IV. MAXIMUM ENTROPY METHOD 

It is an ill-posed problem to extract the spectral func- 
tion from the (imaginary-time) quark propagator. Actu- 
ally, there is an infinite set of spectral functions which 
can reproduce a given correlation function with tolera- 
ble errors. The MEM (l6Hl9Jl considers the probability 
distribution of spectral functions to produce the most 
probable one. The theoretical basis is Bayes' probability 
theorem. The conditional probability of having the spec- 
tral function p(uS) given the correlation function D(t) 
reads 



P[p\DM] 



P[D\pM]P[p\M] 

P[D\M] 



(21) 



where M summarizes all definitions and prior knowledge 
of the spectral function, P[D\pM] and P[p|M] are called 
the likelihood function and the prior probability, respec- 
tively. Since P[D\M] is independent of p(u>), it can be 
treated as a normalization constant. 

According to the central-limit theorem the data D(t) 
are expected to obey a Gaussian distribution: 



P[D\ P M] 



1 



-L[p] 



with 
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(22) 



(23) 



where Zl is a normalization constant, D[p] denotes the 
correlation function reproduced by Eqs. (|16|l or (|17|) given 
the spectral function p(u>), and ^(r) is the variance of the 
error. Maximizing the likelihood function is equivalent to 
X -fitting. 

The construction of the prior probability P[p|M] is the 
central idea of MEM, which expresses the prior in terms 
of the spectral entropy as 



P[p\M(a)} = -L e QS ["' m ] 
Zs 



(24) 



where Zs is a normalization constant and a is an un- 
determined positive scale factor. The Shannon- Jaynes 
entropy S is defined as 



S[p, m] 



dui 



p(ui) — m{ui) — p(u)hi 



m(w) 



,(25) 



where m(tu) is the "default model" of the spectral func- 
tion. Its typical form is a uniform distribution without a 
priori structure assumption [27J, i.e., 



m(w) = m 6(A 2 - us 2 ) 



(26) 



Note that a reliable output from the MEM should be 
insensitive to m and A. If the spectral function is not 
positive definite, one can decompose it in terms of two 
positive definite components, i.e., 

p(«)=p+(«)-p_(w) ) (27) 

Correspondingly, the total entropy is expressed as [29|, |3CJ 

S[p,m] = S[p + ,?7i + ] + S[p-,m-] , (28) 

where m± denotes the default models of p±, respectively. 
Finally, one obtains the total probability distribution 



P[p\DM(a)] (x e 



aS[p,m]—L[p] 



(29) 



The most probable spectral function p a {to) for fixed a can 
be obtained by maximizing P[p\DAI(a)], where usually 
the standard singular-value decomposition algorithm of 
Bryan |16| is adopted. To deal with the scale factor a, we 



follow Bryan's Method [16J. The MEM spectral function 
is defined as 



= / da J Vpp(Lj)P\p\DM(a)]P[a\DAl] 



dap a (u)P[a\DM], 



(30) 



where it is assumed that P[p\DM(a)] is sharply peaked 
around p a (oj), so that the functional integral over p can 
be approximated. In this way, the MEM spectral func- 
tion becomes an average of the p a (u>)'s with respect to a. 
The conditional probability P[a|Z)M] can be evaluated 
using Bayes' theorem as 



P[a\DM]= VpP[p\DM(a)]P[a\M] (31) 

<xP[a\M] fvpe aS[p - m] - L[p] . (32) 

Using the saddle-point approximation and the Laplace 
rule (P[a|M] = const), one obtains 



P[a\DM] oc cxp 



\ k 



a + Xk 



aS[p a ,m] - L[p a ] 



where the Xk are eigenvalues of the following real sym- 
metric matrix in functional space 



A-y = VA 



d 2 L 
dpidpj 



i Pi 



(33) 



Normalizing P[a\DM) and using Eq. ([3U| one finally ob- 
tains pmem- 



V. NUMERICAL RESULTS 

At zero temperature, T = 0, the largest contribu- 
tion to the constituent quark mass comes from DCSB 
which dominates low-energy hadron physics. With in- 
creasing temperature, T > 0, the dynamical quark mass 
decreases, which indicates a partial restoration of chiral 
symmetry. In the chiral limit, there exists a critical tem- 
perature T c where the dynamical quark mass drops to 
zero and chiral symmetry is completely restored through 
a second-order phase transition [3l| . Because of nonzero 
current quark masses, chiral symmetry is not exact. In- 
stead of a second-order phase transition, a crossover hap- 
pens at some pseudo-critical temperature T c which is 
defined by the steepest-descent point for the dynamical 
quark mass. For these two cases, the behavior of the dy- 
namical quark masses with temperature is illustrated in 
the upper panel of Fig. [TJ The (pseudo-)critical temper- 
atures have been indicated as vertical dashed lines. 

Using our model parameters which are able to pro- 
vide a uniformly good description of vacuum properties 
of pseudoscalar and vector mesons with masses < 1 GeV, 
we calculate the dependence of the (pscudo-)critical tem- 
peratures on the interaction width <r, which is shown in 
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FIG. 1. (color online) Upper panel: behavior of the dynam- 
ical quark mass with temperature (for a — 0.5 GeV). The 
black line is the chiral limit, the red line is for the physical 
value of the current quark mass. The black dashed line de- 
notes the critical temperature T c of the second-order phase 
transition in the chiral limit, the red dashed line denotes 
the steepest-descent point for the dynamical quark mass, i.e., 
the pseudo-critical temperature for the physical current quark 
mass. Lower panel: Dependence of the (pseudo-) critical tem- 
perature T c on the interaction width a in our model. 



the lower panel of Fig.Q] T c monotonically decreases with 
increasing a. This behavior is consistent with results ob- 
tained in Ref . [3l| . Remarkably, the critical temperature 
range overlaps well with that obtained by lattice QCD, 
i.e., T c e [0.146, 0.170] GeV [H]. 

Above the critical temperature T c , the quark spectral 
function has been studied by both pcrturbative and non- 
perturbative approaches. At T > 3T C where perturba- 
tion theory (hard-thermal-loop resummation) works, the 
properties of the QGP are dominated by two collective 
excitations: thermal and plasmino excitations [33j. At 
T > T c , experimental obscrvables indicate that nuclear 
matter is a strongly-coupled QGP (sQGP) [34j. In this 
temperature region, perturbation theory fails while the 
nonperturbative DSE approach predicts a novel zero ex- 
citation mode in addition to the normal thermal and plas- 
mino ones |27j . The typical behavior of the quark spec- 
tral function is illustrated in Fig. [2] Here, the spectral 
function is positive definite and each peak corresponds 
to an excitation mode. 
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FIG. 2. (color online) Typical behavior of the spectral func- 
tion at T > T c (following Ref. H3]). 



FIG. 4. (color online) The dynamical quark mass Mo and the 
deconfinement order parameter L p as a function of tempera- 
ture (it = 0.5 GeV, chiral limit). 




-1 
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FIG. 3. (color online) The behavior of the quark spectral 
function at T = 0.8T C (a = 0.5 GeV, chiral limit). 



Below the critical temperature T c , the system is non- 
perturbative because of DCSB and/or confinement. Nev- 
ertheless, the spectral function computed from the solu- 
tion of the truncated gap equation can provide some non- 
trivial information about the system. We first calculate 
the quark spectral function at T = 0.8T C with the inter- 
action width a = 0.5 GeV and in the chiral limit, which 
is plotted in Fig. [3] It is found that the quark spectral 
function exhibits some negative peaks and thus obviously 
L p > 0. Although the physical meaning of those nega- 
tive peaks is unclear, it still makes sense to analyze how 
their behavior changes with temperature. We found that 
the structure of the nonpositive spectral function remains 
unchanged while the residues of the negative peaks, i.e., 
L p , decrease with increasing temperature. Notably, there 
exists a critical temperature Td where L p drops to zero, 
which signals the positivity restoration of the spectral 
function and deconfinement. The calculated behavior of 



Lp is shown in Fig. U in comparison with that of the dy- 
namical mass Mo, which indicates that Td < T c . Next, we 
calculate the dependence of Td on the interaction width 
a both in the chiral limit and with a physical current 
quark mass, which is shown in Fig. [SJ Td monotonically 
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FIG. 5. (color online) The (pseudo-)critical temperature T c 
(black) and the deconfinement temperature Td (red) vs. the 
interaction width a in the chiral limit (dashed lines) and for 
a physical current quark mass (full lines) . 



decreases with increasing ct, and Td is slightly smaller 
than T c . The difference between Td and T c for a physical 
current quark mass is smaller than that obtained in the 
chiral limit. Specifically, when a € [0.4, 0.6] GeV, we have 
T d ~ 94%T C in the chiral limit and T d ~ 96%T C with a 
physical light quark mass. The numerical results are pre- 
sented in Tabled) Our results are consistent with Ref. [18[ 
which also found positivity violations of the Schwinger 
function below T r . 



TABLE I. Critical temperatures of chiral symmetry restora- 
tion T c and deconfinement Td for different parameters (di- 
mensionful quantities reported in GeV, A = T c — Td). 
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rpO 
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A m/ T m 


0.4 


0.165 


0.155 


6.1% 


0.168 


0.160 


4.8% 


0.5 


0.142 


0.131 


7.7% 


0.144 


0.136 


5.6% 


0.6 


0.126 


0.122 


3.2% 


0.128 


0.127 


0.8% 


avg. 


0.144 


0.136 


5.7% 


0.147 


0.141 


3.7% 



By defining a confinement scale r a = 1/er, it is appar- 
ent that both T c and X^ increase with increasing r a , or 
T c ,d oc r CT . Considering that the difference between T c 
and Td is just several MeVs, while the systematic un- 
certainty introduced by our approximations is certainly 
larger, it is reasonable to claim that chiral symmetry 
restoration and deconfincment coincide at nonzero tem- 
perature and zero chemical potential. 



VI. SUMMARY AND REMARKS 

At nonzero temperature and zero chemical potential, 
we computed the quark spectral function via the MEM 
from a solution of the QCD gap equation. For the lat- 
ter, we used a rainbow interaction kernel which phe- 
nomenologically models recent results from DSE and lat- 
tice QCD. As a criterion for the positivity violation and 
restoration of the quark spectral function, we proposed 
an order parameter L p which is directly related to the 
integral of the spectral function's negative part and ob- 
viously vanishes for positive definite spectral functions. 
We indeed found that L p > at low temperature while 
L p = at high temperature, i.e., there exists a criti- 
cal temperature Td where L p drops to zero. Here, the 



positivity of the quark spectral function is restored and 
quarks become asymptotic particles for T > Td- There- 
fore, we conjecture that the deconfinement phase transi- 
tion happens at T = Td- Using our model setup, which 
can uniformly well describe vacuum properties of pseu- 
doscalar and vector mesons with masses < 1 GeV, the 
critical temperature of deconfincment comes out slightly 
smaller than that of chiral symmetry restoration, i.e., 
Td ;$ T c . Within the systematic uncertainties of our ap- 
proach, however, it is reasonable to conclude that chiral 
symmetry restoration and deconfincment coincide. 

It is generally expected that nuclear matter has a rich 
phase structure because of the interplay between DCSB 
and confinement at nonzero chemical potential. It would 
therefore be interesting to extend the present study to 
the case of nonzero chemical potential. The difference 
between T c and Td could then be more pronounced and 
give rise to the so-called quarkyonic phase [35], [36| . 
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